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We consider analytically and numerically an anisotropic spin-| delta-chain (saw¬ 
tooth chain) with nearest-neighbor ferromagnetic and next-nearest-neighbor antifer¬ 
romagnetic interactions. For certain values of the interactions a lowest one-particle 
band becomes flat and there is a class of localized-magnon eigenstates which form 
a ground state with a macroscopic degeneracy. In this case the model depends on 
a single parameter which can be chosen as the anisotropy of the exchange interac¬ 
tions. When this parameter changes from zero to infinity the model interpolates 
between the one-dimensional isotropic ferromagnet and the frustrated Ising model 
on the delta-chain. It is shown that the low-temperature thermodynamic properties 
in these limiting cases are governed by the specific structure of the excitation spec¬ 
trum. In particular, the specific heat has one or infinite number of low-temperature 
maxima for the small or the large anisotropy parameter, correspondingly. Numeri¬ 
cal calculations of finite chains demonstrate that this behavior is generic for definite 
values of the anisotropy parameter. 



I. INTRODUCTION 


Quantum many-body systems with a single-particle flat-band have attracted much at¬ 
tention [l-7|. Frustrated quantum spin systems represent examples where flat-band physics 
may lead to new interesting phenomena such as a nonzero residual ground state entropy. 


extra low-temperature peak in the specihc heat etc j7-lll|. An interesting and typical ex¬ 


ample of such flat-band system is the s = ^ delta or sawtooth Heisenberg model consisting 
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FIG. 1: The A-chain model. 


of a linear chain of triangles as shown in FigHJ The interaction Ji acts between the apical 
(even) and the basal (odd) spins, while J2 is the interaction between the neighboring basal 
spins. The Hamiltonian of this model has a form 

N 


^ / 1 \ 
+ (^S2j_iS2j+i + s|j_iS2i+l + ^2{s2i-lS2i+i - -)j 


( 1 ) 


where are s = | operators, Ai and A2 are parameters representing the anisotropy of 
basal-apical and basal-basal exchange interactions respectively, N is the number of sites. 
For the periodic boundary conditions (PBC) si = sat+i. The constants in Eq.([l]) are chosen 
so that the energy of the ferromagnetic state with is zero. 

The ground state of the isotropic Heisenberg model (Ai = A2 = 1 ) with both antifer¬ 
romagnetic interactions (Ji, J2 > 0 ) (AF-delta-chain) has been studied as a function of 


parameter 


J2 

Ji 


12 


14 l |. Remarkably, for a special choice of the interaction values Ji = 2J2 


the lower one-magnon band is dispersionless and the excitations in this band are localized 
states. The localized nature of the one- magnon states is a base for the construction of 
multi-magnon states because a state consisting of k independent (non-overlapping) local¬ 
ized magnons is an exact eigenstate. Such construction is possible up to A; < ^ and these 
states form the ground state manifold at the saturation magnetic held. The ground state 
and low-temperature properties of the AF delta-chain with Ji = 2 J2 have been studied in 


detail in Refs. jr-llOj]. Typical features related to the localized magnon states are the zero- 
temperature magnetization plateau and the magnetization jump, residual entropy and extra 
low-temperature peak in the specihc heat. 

In contrast to the AF delta-chain the same model with Ji < 0 and J2 > 0 (F-AF 
delta-chain) is less studied though this model is used as a model of several compounds 
such as [Cu{bpy)H20] x [Cu{bpy){mal)H20]{Cl04)2 containing magnetic ions Q- 
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. Similarly to the AF delta-chain in this model the localized states exist if |Ji| = 2J2. 


It is known 


18l | that the ground state of the F-AF isotropic delta chain is ferromagnetic 


for Oi = < |. In Ref.[l 8 | it was argued that the ground state for a > ^ is a special 


ferrimagnetic state. The critical point a = ^ is the transition point between these two 
ground state phases. The isotropic F-AF delta-chain at the transition point ol = \ has been 


studied in Ref. 


19 | . It was shown 19 | that in addition to the multi-magnon conhgurations 


consisting of isolated magnons the special states with overlapping magnons (localized multi- 
magnon complex) exist and all of them are exact ground states at zero magnetic held. So, the 
ground state degeneracy in this model is even larger than for the AF delta-chain. Another 
difference between two isotropic models concerns the energy gaps between the ground state 
and the excited states. In the AF delta chain these gaps are hnite while in the F-AF 
model the gaps for /c-magnon states decrease rapidly with the increase of k. As a result the 
contribution of the excited states to the thermodynamics can not be neglected. 

It is interesting to study the inhuence of the anisotropy of exchange interactions in the F- 
AF delta-chain on the ground state properties and on the low-temperature thermodynamics. 
As it will be shown there is a special line in (Ai,A2) plane on which the localized magnons 
are exact ground states in zero magnetic held. The main aim of this paper is to study the 
F-AF delta-chain on this line. We will demonstrate that the behavior of the model on this 
line has non-trivial peculiarities 

The paper is organized as follows. In Section II we derive the conditions on model pa¬ 
rameters that provides localized magnon eigenstates and, therefore, macroscopic degeneracy 
of the ground state, which is calculated in Section III. In Section IV we study the low- 
temperature thermodynamic of the system both analytically and numerically. In Section V 
we give a summary of our results. 


II. ONE-MAGNON STATES 


We begin the study of the anisotropic F-AF chain described by Eq.([T]) with the one- 
magnon spectrum over ferromagnetic state. Two branches of states with — 1 are 

given by 


E^{q) = Ai - 


1 


a(A2 — cos q) ± \Ja^{^2 — cosqY -|- 2(1 -|- cosg) 


( 2 ) 
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here and further we put Ji = — 1 . 

At a dehnite value of a = ao with: 


ao — 


^2(1 + A 2 ) 

the lower band becomes dispersionless with the energy 

1 


(3) 


e: = Ai 


2Q;n 


(4) 


We note that the value of Oq does not depend on Ai but the energy e: does. The disper¬ 
sionless one-magnon states correspond to localized states which can be chosen as 


\F) = {S2i + —S2i+1 + S2i+2) \F) i = (5) 

ao 

where \F) is the ferromagnetic state with all spins up, is on-site spin lowering operator 
and n = N/2. 

The wave function ipi |F) is localized in a valley between Tth and {i + l)-th triangles. 
The wave functions Eq.(l5]) are eigenfunctions of at a = ao- To prove it let us represent 
the Hamiltonian H as a sum of local Hamiltonians 


= ( 6 ) 

i 

where Hi is the Hamiltonian of the Tth triangle which is 

Hi = - ^ (s2i+5S2i + ^2i+5^2i + ^l{s2i+5^2i ~ 7^) 

< 5=±1 ^ ^ ^ 

+a (^S2i_iS2i+i + sli-isli+l + ^2{s2i-lS2i+i - -)^ (7) 

It is easy to check that at a = cto 

{H, + H,+,)if,\F)=e0,\F) ( 8 ) 

and Hj(pi |F) = 0 for j 7 ^ i, i - 1 - 1 . 

Therefore, 

Hip, \F) = ep, \F) (9) 

Because the one-magnon wave function Eq.(|S]) is localized it is possible to construct the 

states with k independent (non-overlapping) localized magnons for k < ^ with the energy 
Ek = ek. 
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Thus we found that model Eq.(IT]) for definite choice of parameter a = Uq given by Eq.(|3]) 
has localized magnon eigenstates Eq.(|5]) with the energy £ Eq.(|l]). If parameters Ai and A 2 
are chosen so that £ < 0 (A 2 > 2A^ — 1), all the states composed oi k < ^ independent 
localized magnons are the lowest ones in the corresponding spin sector — k. It turns 

out that the lowest state in the case e < 0 lies in the sector = 0 |^. The magnetic 
properties of the F-AF delta-chain at a = ao and e < 0 are similar to those for the AF 
isotropic delta-chain with A = 2 J 2 . In particular, the ground state magnetization curve has 
the plateau and the jump, and all the states composed of independent localized magnons 
form the macroscopically degenerated ground state at the saturation magnetic held. 

If £ > 0 (A 2 < 2A‘l — 1) the localized magnons are exact eigenstates as well, but they 
are not the lowest eigenstates. Moreover, if e: > 0 the ground state is ferromagnetic with 
Sfatal — prove it let us consider eigenvalues of the triangle Hamiltonian with a = uq. 

The spectrum of each local Hi consists of four levels, all of them are two-fold degenerated 
over Sz'- 

El 

E 2 
E2,,i 

It follows from Eq. flT0|) that all eigenvalues for Sz = in Eq. flT0|) are positive at e > 0. 
Since the local Hamiltonians of neighbor triangles do not commute with each other the 
ground state energy Eq of H satishes an inequality 

-Fq > ^ Eoi = 0 (11) 

i 

where Em = 0 is the ground state energy of i-th. triangle with Sz = ±|. The inequality flTT|) 
turns in an equality only if all triangles have 5*^ = | or = — | simultaneously. Then, the 
states with Sfatal = two ground states only. 

Thus, we conclude that the ground state with Sfatai = 0 is realized if £ < 0 while for e > 0 
the ground state is ferromagnetic. However, if e = 0 the lowest eigenvalues with 5*2 = ±| of 
the local Hamiltonian becomes zero and it indicates the possibility of a signihcant increase 
of the ground state degeneracy. 

The condition e = 0 means that model parameters Ai and A 2 are coupled and we 


= 0 , 52 = ± 


e 

2 ’ 

1 


= ± 


2 '““ + !+2a„ 



(o^O + - + 7^-— 2 q!o£, 

Z ZOlQ 


s, = ±- 


( 10 ) 
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parameterize them by the anisotropy of the basal-apical interaction Ai 


Aa = 2Al-l 


( 12 ) 


Hamiltonian ([T]) of the anisotropic F-AF delta-chain with this choice and at a = Oq has 
a form 


» = - E («+i + 4»r+, + A, - i)) 

+^E + 4-.4.+1 + (2a; - 1) (4-.4i+i - J 


(13) 


Model ffT3|) is the main object of the following study. It has macroscopic degeneracy 
of the ground state and nontrivial thermodynamic properties. When the parameter Ai is 
changed from Ai = 0 to Ai = oo the model interpolates between the isotropic ferromagnetic 
chain on the basal sites and the Ising model with equal but opposite in sign basal-apical 
and basal-basal interactions. In the case Ai = 1 model (IT^ reduces to the isotropic F-AF 
delta-chain studied in R ef.IlQ l. 

n rr 

As noted in Ref. |18l. Il9| . Hamiltonian (IT^ with Ai = 1 describes the model at the 
transition point between the ferromagnetic and the ferrimagnetic ground states. In this 
respect model flT^ as a function of Ai describes the transition line of more general model, 
where the parameter a is not put as a = 2 ^- This model reads: 

1 ' 


^ = -E 


r.a; X 


-I 


+ ^1 1 *^*i+l ~ 


+“E (^' 522 - 1 * 22+1 + *22-1*Ei + (2^1 - 1) (^*22-1*22+1 “4)) 


(14) 


The ground state phase diagram of model flTT)) in (Ai,q;) plane obtained by numerical 
calculations is shown in Figj2j The curve a = 2 ^ corresponds to the transition line between 
different phases shown in FiglJl We note that the study of the behavior of the model in 
these phases is out of scope of the present paper and will be given elsewhere. Here we focus 
on Hamiltonian (1X31) . 


III. THE GROUND STATE DEGENERACY 


In this Section we study the ground state of Hamiltonian (lT3il . 
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FIG. 2; The ground state phase diagram of model ([HI). 


As follows from Eq. llTTl) the ground state energy of this Hamiltonian is zero. There are 
n one-magnon states Eq.(|5]) in the spin sector 5^^^; = y — 1. These states form a complete 
nonorthogonal basis, have zero energy and, therefore, belong to the ground state manifold. 

It is evident that the pairs of the isolated magnons (|^ — j | > 1) are the eigenfunc¬ 
tions with zero energy of each local Hamiltonian (and, therefore, of total Hamiltonian) in 
the spin sector = y — 2. Similarly, eigenfunction composed of k isolated magnons is 
the exact eigenfunction of the ground state in the sector = y — fc. However, such set 


of functions contains not a 
As it was shown in Ref. 


1 of the states with zero energy in fc-magnon sector with k > 1. 


19| for the isotropic F-AF delta-chain (Ai = 1) there are also the 


ground state eigenfunctions of another type, and this type of functions holds in model flT^ 
with 0 < Ai < cxo as well. As an example, we consider two-magnon eigenstates. For k = 2 
along with the pair of isolated magnons we can write the exact two-magnon state as 


+ Bipi -F ^i+i) \F) 


(15) 


where B = 2A^ — 1. 

It is easy to check that function flT^ is the exact eigenfunction with zero energy for the 
local Hamiltonians ifj, Hi+i and and for the others ones. For Ai = 1 this function 


reduces to that considered in Ref.[l9|]. We note that the function flT^ contains overlapping 
magnons. The complete manifold of ground states in the sector = y — 2 consists of 
pairs (n = A/2 is assumed to be even) of independent magnons and n eigenfunctions 










ffT^ . Thus, the ground state degeneracy in this spin sector is 

The construction of the wave functions of type Eq. ffTSj) can be extended for k > 2. 
According to the results of Ref. 19(] the total number of the ground states of the isotropic 
F-AF chain for hxed = {n — k) {0 < k < n) is 


GJk) = Cl 

Gn{k) = + h.n, 


, n 
0<k<- 


n 

— < k < n 
2 - - 


( 16 ) 


The total degeneracy of the ground state W is 


IF = 2" + + 1 


(17) 


so that W = in the limit iV —)■ oo. 

V 'TT 

The analysis similar to that for the isotropic model can be carried out for the anisotropic 
delta-chain Fq.([T3]). It turns out that the degeneracy of the ground state given by Fqs. (lT6|l - 
cm is valid for all 0 < Ai < cxd except two special cases Ai = l/\/2 and Ai = 1/2 which 
are considered below. 

In the case Ai = l/\/2 the anisotropy A 2 vanishes and the model reduces to 




i '’i+1 




y 

2i+l 


(18) 


y /2 * ^/2 

In this special case the states containing neighbor localized magnons like (/’jt/’i+i are exact 
ground states (R = 0 in Fq. flTHl) ). which is not valid in general case. We note that similar 
type of exact eigenstates arises in the Hubbard model on the delta-chain in which neighboring 
valleys can be occupied by the localized electrons with identical spins j^. Q . The presence 
of such exact states for model Fq. dTSD leads to the increase of the ground state degeneracy, 
which is 


W = 2 *^ - + 1 


n 




(19) 


In the case Ai = 1/2 after rotation in the XY plane —>■ (—1)^3^^ the Hamiltonian 

takes the form 


^ = E ( 


E 4-i4+i + 4-i4+i-o4-h 


2i+l 


( 20 ) 


That is basal-apical and basal-basal interactions are exactly equal in this case. This 
fact raises symmetry of the system: isosceles triangles in general case becomes equilateral 
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triangles in this special case. The increased symmetry results in the additional degeneracy 
of the ground state, which is 

n/ = 2" (^ + l) + + 1, Ai = i (21) 

In the limit Ai = 0 the model reduces to the quantum ferromagnet on basal spins and 
independent n apical spins. Therefore, the ground state degeneracy is 


IT = (n + l)2”, Ai = 0 


( 22 ) 


where the factor (n + 1) comes from the degeneracy of the ferromagnetic state over total S^. 

The special case Ai —)■ cx) will be studied in detail in Sec.IV. Here we give only the results 
for the ground state degeneracy: W = 3” + 1. 

It is interesting to note that the ground state degeneracy of the anisotropic delta-chain 
with the open boundary conditions (OBC) is the same as for the isotropic model for all 
value of 0 < Ai < CX). According to the results of Ref. the ground state degeneracy of 
the open chains with add N sites is 

N + l 


WoBC = ‘2,’^~\n + l), 


n = 


(23) 


All the above presented expressions for the ground state degeneracy have been conhrmed 
by exact diagonalization (ED) calculations of hnite chains. 

The exponential degeneracy of the ground state results in the residual entropy Sq = 
\n{W)/N. Though the numbers of degenerated states in the general case Eq. flTT)) and in 
special cases Eqs. (lT^ . fl^ .fl2 ^ . fl2^ are different, they yield the same result for the residual 
entropy in the thermodynamic limit N ^ oo: 

So = ^ In 2 (24) 


The difference in the numbers of degenerated ground states reveals itself in the corrections 
~ \n{N)/N, that vanishes in the thermodynamic limit. 


IV. LOW-TEMPERATURE THERMODYNAMICS 
A. Ising model 

In this Section we study the low-temperature behavior of the anisotropic F-AF delta- 
chain on the transition line. We begin with a limit Ai —)■ cx) when model flT^ reduces to 
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the Ising model, the Hamiltonian of which has a form: 


H, = 


E 


(s. 


i^i+l 


- j) + E(4-i4+i- 


(25) 


where we introduced the dimensionless magnetic held h. 

The partition function of this model can be obtained using a transfer-matrix method and 
given by 

Z = XI+ \l (26) 


where eigenvalues of the transfer-matrix are 

Ai = 1 + 2 cosh + exp 

A 2 = -1 + exp (27) 

The ground state of model fl2^ at h = 0 has zero energy and the ground state degeneracy 

Gnik) in the spin sector = {n — k) {Q < k < n) can be found as coefficients in the 

expansion of Z in powers of exp(—h/T). As a result Gn{k) is given by 

Gn(k) = ■£ + h.n ( 28 ) 

m=0 

The total ground state degeneracy is 

n-l 

IT = 2 ^ Gnik) + Gnin) = 3^ + 1 (29) 

fc =0 

The residual entropy per site Sq of the Ising model fl2^ at —)■ cx) equals 

So = ^ ln3 (30) 


Using Eqs. (l2^ and (1271) we can obtain all thermodynamic quantities (in the thermody¬ 
namic limit only largest eigenvalue Ai survives). In particular, the specihc heat per site is 
given by 


G = 


2T2 


3 exp + exp(-^) 


(31) 


The specihc heat as a function of temperature has a typical broad maximum around 
T ~ 0.5 and exponential decay for T —)■ 0. The zero-held susceptibility per site x is given 


by 


X = 


T 


3 + exp (-y:) 


(32) 
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and it behaves as x ~ 3 ^ for T —)■ 0. 

Now let us consider the generalization of the Ising model fl2^ where basal-apical and 
basal-basal interactions are different: 


Hj = 


El 


SiS. 


i+1 


7) + (1 + 7) E('® 2 i-l' 52 i+l - 7) 




(33) 


The ground state of model fl33p is ferromagnetic for 7 < 0 and the ‘antiferromagnetic’ 
with = 0 on the basal subsystem for 7 > 0. The Ising model (125|1 (7 = 0) describes 
the transition point between these phases. Generally speaking, the Ising model fl3^ with 
7 7 ^ 0 is not any limiting case of the initial model flT3|) . Nevertheless, it is useful to study 
model fl33|l with 7 > 0 because on the one hand it has exact solution and on the other hand 
some properties of the thermodynamics in the AF phase, especially in the vicinity of the 
transition point, are inherent in model flT^ as well. 

The eigenvalues of the transfer-matrix are 


Ai 2 = cosh 


T 


+ exp ( -- ) ± 


\ 


1 - 1 - cosh 


T 


2 exp ( ^ +cosh 


- 1 


(34) 


The partition function at S> 1 is 


Z = 2^ exp 


7n\ 

2t) 




2-^7 

2T 


(35) 


The ground state degeneracy W = 3" for 7 = 0 is partially lifted up to 2" for 7 > 0. The 
degeneracy 2 ” is related to an independence of the ground state energy for 7 > 0 on the 
spin conhguration of the apical subsystem. 

The temperature dependence of the specihc heat obtained from Eq. (|35il has a form 


(36) 


The specihc heat as a function of T is shown in Figj3]for 7 = 0.001. In comparison with 
the case 7 = 0 the dependence C{T) has an additional low-temperature maximum at T — 
Such two-peak form of the temperature dependence of the specihc heat can be explained as 
follows. The spectrum of the Ising model (13^ with 0 < 7 <C 1 has two energy scales. One 
of them is set by (4"" — 3”) states with energies E ~ 1 and another one is set by (3*^ — 2 ”) 
low-energy states that split-oh from the ground state at small 7 . According to Eq. fl3^ 
the density of the low-energy states having the energy E = k'y/2 is p{k) = with a 
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T 

FIG. 3: The specific heat and entropy as the function of T of model (1331) for 7 = 0.001. 

maximum aX k = ^. This part of the energy spectrum is responsible for the appearance 
of the low-temperature maximum in C{T). The two-scale form of the low-energy spectrum 
causes the peculiarities of the other thermodynamic quantities. For example, the entropy 
per site s(T) has stair-step-like temperature dependence (see Figl3]). In the regions T -C 7 , 
7 -CT-CI, 1<^T the entropy per site s{T) behaves as s ~ 1 In 2, s ~ | ln3, s ~ | ln4, 
correspondingly. The quantity yT has similar stair-step-likes behavior. 

B. F-AF delta-chain with Ai 1 

When the parameter Ai is large it is convenient to normalize Hamiltonian ffT3|) as H/Ai 
and to write it in the form 

-^H = Hj + V, + V2 (37) 

where Hj is the Ising Hamiltonian fl2^ at h = 0 and Vi and V 2 have the forms 

V 2 = 2g‘^ - ^2i-1^2i+l + (38) 

with small parameter g = l/( 2 Ai). 

At 5 ^ = 0 the ground state of Hamiltonian fl3711 is 3”'-fold degenerated and there are 
(d"" — 3") states with E >1. The terms Vi and V 2 lift the degeneracy for each spin sector, 
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but partly only: the part of the ground state levels remains degenerated with zero energy 
(the number of them is given by Eq. ffT7|) i and other ones move up. As a result the ground 
state degeneracy of model flTTl) is IV ~ 2 ” in comparison with Wj = 3” at 5 ^ = 0 . At 5 ^ -C 1 
the levels of Hamiltonian (l37|l which split off from the ground state form a set of low-lying 
excitations determining the low-temperature thermodynamics. 

The calculation of the spectrum of these states is very complicated problem and we begin 
with the one-magnon subspace = y — 1. In the one-magnon subspace the number 

of degenerate ground states of model fimi is n for both cases g = Q and 7 ^ 0 , i.e. the 
perturbation V = V + V does not lift the degeneracy. In the two-magnon sector the ground 
state degeneracies of Hi and of H are C'^+i and correspondingly. It means that n states 
split off from the ground state. In fact, these states are two-magnon bound complexes and 
their energy found in the second order m g is E = g'^. 

The calculation of the PT in g for A: = 3 is rather cumbersome and we give the hnal result 
only. The set of (C^ - 1 - n{n — 1)) three-magnon states which are degenerated at 5 ^ = 0 splits 
into three subsets. One of them contains ground states with E = 0. There are n{n — 2 ) 
states with E ^ g"^. These states consist of the two-magnon bound complex and the isolated 
localized magnon. The third subset represents n three-magnon bound complexes with the 
energy E = 3g‘^. 

The analytical computation of the spectrum of low-lying excitations for fc > 3 is com¬ 
plicated problem and our further conclusions are based on numerics. According to nu¬ 
merical data the structure of the spectrum of low-lying fc-magnon states is similar to that 
for the case k = 3. The perturbation V = Vi + V 2 splits the ground state manifold into 
fc-subsets: the ground states with the energies E = 0] A:-magnon bound complexes with 
E ~ ^ 2 (A:- 1 ). the states consisting of (k — l)-magnon bound complex and one isolated 
magnon (E ~ the states consisting of {k — 2)-magnon bound complex and two 

isolated magnons {E ~ and so on. The highest subset of excitations has the ener¬ 

gies E ~ g^. Thus, the low-lying excitations in the sector with = y — A: is distinctly 
divided into the parts with the energies E ^ g'^, E ^ g^,... E ^ 

Taking into account all states with all possible values of we found that the total 

spectrum of model fl37D can be rank-ordered in powers of small parameter g^ and it has a 
multi-scale structure. The ground state degeneracy is given by Eq. ffT7|l . The lowest excited 
states for hnite chain composed of n triangles has the energy E ~ ^ 2 n -2 j^eans that 


14 



E 

FIG. 4: The integrated density of states of normalized Hamiltonian (1371) for g = Q.l and N = 16. 
the gap is exponentially small. 

The distribution of the energy levels for = 16 and = 0.1 is shown in FigJU As it is 
seen in FigJU the spectrum is distinctly divided into the parts. Each part of the spectrum 
behaves as E ~ as written in FigJU This fact was conhrmed numerically by comparison 
of the energies for g = 0.1 and g = 0.125. We note that for = 16 and g = 0.1 the 
value of the smallest excitation E ^ becomes very small {E 10 — 10 and 

it is indistinguishable from the ground state, so only the excitations with E > 10“^^ are 
accessible by the ED calculations and represented in FigJU 

Such structure of the spectrum determines the characteristic features of the low- 
temperature thermodynamics. To study the thermodynamics of model flT3|) we use the 
exact diagonalization (ED) of hnite delta-chains with PBC up to A^ = 20. In FigsJS] and J7| 
we represent the data for the specihc heat C{T) and the entropy s{T) (both per site) for 
Ai = 5 {g = 0.1) obtained by the ED for A^ = 16 and N = 20. The temperature depen¬ 
dence of the specihc heat shown in Fig J5] exhibits numerous maxima the presence of which 
can be explained as follows. As it was shown in Sec.lVa the dependence C{T) of the Ising 
model has one low-temperature maximum which is related to the part of the spectrum with 
the energy E ~ 7 ^. Similarly, multi-scale structure of the spectrum of Eq. fJTTJ) for small 
g leads to the dependence C{T) with many maxima related to the corresponding parts of 
the spectrum. This is conhrmed by the following calculation. Let us select the part of the 
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T 

FIG. 5: The dependence of specific heat on the temperature of normalized Hamiltonian (I37p for 
g = 0.1 and N = 16,20. 
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FIG. 6: The dependence of specific heat on the temperature of normalized Hamiltonian (j37p for 
g = 0.1 and N = 6,8. 

spectrum with E g^, remove the high-energy part with i? ~ 1, put the energy of all 
remaining low-lying states to zero (as the ground state) and calculate the contribution of 
such deformed spectrum to the specihc heat. It turns out that this contribution perfectly 
describes the hrst low-temperature peak in C{T). Similar procedure for the parts of the 
spectrum with E g^, E and so on reproduces very well corresponding peaks in C{T). 
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The chain consisting of n triangles has n peaks in the dependence C{T). Thus, the 
number of peaks grows linearly with the system size N. It is illustrated in Figj 6 l where the 
specihc heat C{T) is shown for N = 6 and N = 8 with g = 0.1. As it can be seen from Figj 6 ] 
the number of peaks is three (four) for = 6 ( 8 ). 

Since the m-th peak arises from the part of the spectrum with E ~ the maximum 
of the m-th peak takes place at 

Tm ~ {cgf^ (39) 


with some constant c. 

If the value g is small the corresponding temperature Tm is very small as well. For this 
reason we could not represent in FigJS] all feasible peaks for = 16 and A^ = 20 because 
the temperatures Tg and T 7 for = 0.1 are out of the accuracy of the ED calculations. 

As it follows from the above the lowest peak should occur at T„ ~ (cg)^'^. This allows us 
to write the hnite-size scaling parameter in the form 


InT 


(40) 


The dependence of the entropy per site on temperature for small g has stair-step behavior 
as shown in FigJTl These stair steps he in between the limiting values s = ln2 for T —>■ 00 
and s = i In 2 for T —)■ 0. 

As it can be seen from Figsj5]and[7]the data of C{T) and s(T) for A^ = 16 and A^ = 20 
deviate from each other for T < Tq = 10“® but they are indistinguishable for T > Tq. 
This indicates that the obtained finite-size data correctly describe the thermodynamic limit 
for temperatures T > Tq. For example, at least three low-T maxima for C{T) and three 
stair-steps in the dependence s{T) remain the shape at A^ —)■ cxd as in FigsJS] and [3 The 
levels of the entropy s{T) on these three steps testify that the number of states in sectors 
is: 4*" states with E ~ 1; S’" states with E ~ g"^] 2.62*" states with E ~ 51 ^; 2.45*" states with 
E ~ g^. 

The deviation of the data for A^ = 16 and A^ = 20 in the region T < Tq means that the 
hnite-size effects become essential for T < Tq and the correct description of the thermody¬ 
namics in this temperature region needs more large systems. Nevertheless, the multi-scale 
structure of the spectrum for Ai ^ 1 will lead certainly to the existence of n maxima in 
C{T) and n stair-steps in s{T) for chains with n triangles. 
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FIG. 7: The dependence of entropy on the temperature of normalized Hamiltonian (|37p for g = 0.1 
and N = 16, 20. 


The temperature dependence of xT {x is the uniform magnetic susceptibility per spin) 
has a stair-like behavior similar to the dependence s(T). The numerical data show that the 
quantity yT tends to a hnite value depending on At at T —)■ 0 (see Figj8]). This hnite value 
is related to an average of square of over the degenerated ground state as 


XT={{SX)^ 

Evaluating this average with use of Eq. flTbj) we obtain 



(41) 


N 

{xT)n = ^, N ^oo (42) 

According to this equation xT is proportional to N for At 3> 1. It means that xT diverges 
at T —)■ 0 in the thermodynamic limit. The behavior of yT for large N and small T depends 
on the scale variable given by Eq. (l40|) . We will determine the asymptotic behavior of yT 
at T —)■ 0 using the following estimations. We assume that the stair-step-like behavior as in 
FigJH] remains in the limit T ^ 0, so that all steps becomes equal in this limit. Then, for the 
determination of the low-temperature limit of the susceptibility we need to find the height 
and width of the steps. According to Eq. fl42l) the total height is (yT)„ = This means 
that each additional triangle in the system leads to the additional stair-step of the height 


hstep = {xT)n+l - {xT)n 


1 


(43) 
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FIG. 8: The dependence of magnetic susceptibility on the temperature of normalized Hamiltonian 
daZI) for 5 = 0.1 and N = 16, 20. 


The width of each stair-step can be estimated as the distance between neighbor peaks in 
the specihc heat, which is 

T 

(44) 


“.lep = >n = 2 In {eg) 

2 m+l 


The low-temperatnre limit of the snsceptibility can be calculated as envelope of all stair¬ 
steps, so that for xT we obtain 

xT = ^^lnT (45) 

^ step 

Substituting Eqs. flT3|l and (1441) into Eq. fl45|) we obtain the low-temperature asymptotic 
of the susceptibility 

X =-- (46) 

^ 48Tln(c^2) 1 ) 

Susceptibility Eg. (1461) can be rewritten in the scaling form which takes into account the 
hnite-size effects: 

JV 

(47) 




Here the scaling function f{y) of the scaling variable Eq. fl40|) has the limits f{y) = y for 
?/ -C 1 and f{oo) = 1. Thus, in the thermodynamic limit x ^ at T —)• 0. 
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C. F-AF delta-chain with Ai ^ 1 


In the limit Ai = 0 model flT^ reduces to quantum ferromagnet on basal sites and non¬ 
interacting apical spins. For Ai -C 1 it is convenient to normalize the Hamiltonian (IT^ in 
a form 

AiH = Ho + V (48) 

where Hq is the Hamiltonian of the isotropic ferromagnetic basal spin chain and V is the 
basal-apical interaction 


Ho 

V 


1 

2 



X X 

2i-1^2i+l 


+ S 


y oV 

2i-l*2i+l 


Szi-lSzi+l 



-A, E (sfsf+i + 44+i) + 


(49) 


The ground state with zero energy of Ho is (n-l-l)-fold degenerated and all eigenfunctions 
of Ho do not depend on the spin state of the apical subsystem. Therefore, the degeneracy 
of the ground state of model fl48|) at Ai = 0 in the spin sector = n — k is 

G„(k) = y C” (50) 

m=0 

and the total number of the ground states at Ai = 0 is 

IT = 2”(n+l) (51) 


The perturbation V lifts the degeneracy partly and the number of levels which split off 
from the ground states for hxed is Y.^2o C*™- These levels form the spectrum of low- 
energy excitations of Hamiltonian fHHjl for Ai -C 1. In the second order in Ai the lowest 
energy (the gap) in the spin sector = n — k for 0 < fc < | is 

E(k) = A;(1 - hhA) (52) 

n 

For example, the gap for the one-magnon excitations is F^(l) = A^. This agrees with the 
energy of the upper branch of the one-magnon states Eq.(l2]) at g = tt. As to the gaps in 
the spin sectors ^ < k < n, the numerical calculations of hnite chains show that they are 
also proportional to Af. Therefore, the spectrum of model flT8|l at Ai -C 1 has two-scale 
structure with the energies E ~ A^ and E ~ 1. It is illustrated in Fig J9] where the integrated 
density of states is shown for A = 16 and Ai = 0.05. The two-scale spectrum for Ai 1 
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FIG. 9: The integrated density of states of normalized Hamiltonian (148 p for Ai = 0.05 and N = 16. 



10-'= 10-5 ig.4 10-3 0.01 0.1 1 10 

T 

FIG. 10: The dependence of specific heat on the temperature of normalized Hamiltonian (|48l) for 
Ai = 0.05 and N = 16,20. 

is in contrast with the multi-scale spectrum for Ai > 1. Such two-scale structure of the 
spectrum leads to the temperature dependence of the specihc heat with two peaks as shown 
in Fig{Tn]for Ai = 0.05. This behavior is in contrast with multi-peaks form of C{T) for 

Ai > 1. 
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FIG. 11: The dependence of specific heat on the temperature of normalized Hamiltonian (1481) for 
Ai = 0.5 and N = 16,20. 

D. F-AF delta-chain with intermediate valne of Ai. 

As follows from the results obtained above, the structure of the spectrum of model flT^ as 
well as the low-temperature thermodynamics is qualitatively different in the limiting cases 
Ai 3> 1 and Ai 1. When the parameter Ai is neither large nor small the situation is less 
clear. 

Let the parameter Ai moves down from the limiting case Ai S> 1. Then the ordered 
partition of the spectrum is smeared and the multi-scale structure of spectrum becomes 
less obvious. Nevertheless for Ai > 1 the low-temperature behavior of the thermodynamic 
quantities remains qualitatively similar to that for Ai S> 1. For example, the temperature 
dependence of C{T) for the F-AF delta-chain with Ai = 2.5 and Ai = 1 is shown in 
Fig{T2]and Fig{T3l The specihc heat is characterized by the existence of the low-T maxima, 
though they are not so distinctive as in the case Ai = 5 (Figj5]). The reason of this fact is 
that the hnite-size effects become more pronounced when Ai decreases. For example, the 
temperature Tq below which the data for N = 16 and N = 20 start to deviate from each 
other is Tq ~ 10“®, 10““^, 10“^ for Ai = 5, Ai = 2.5, Ai = 1, correspondingly. Nevertheless, 
we expect that the number of low-T peaks of C{T) is proportional to N. The behavior 
of other thermodynamic quantities such as the entropy s(T) and xT(T) for Ai = 2.5 and 
Ai = 1 is qualitatively similar to those for Ai S> 1. This suggests that model flT^ with 



22 



IQ-a j^q-7 iq -6 20-5 lO* 10'^ 0.01 0.1 1 10 

T 

FIG. 12: The dependence of specific heat on the temperature of normalized Hamiltonian (1371) for 
g = 0.2 (Ai = 2.5) and N = 16, 20. 



FIG. 13: The dependence of specific heat on the temperature for Ai = 1 and N = 16,20. 

Ai ;:§> 1 is the generic one for all values of the anisotropy parameter in the range Ai > 1. 

However, for 0 < Ai < 1 the low-T behavior of the thermodynamics changes qualitatively. 
As it can be seen from FigHUJand FigJTT]the specific heat C{T) has more complicated form. 
In addition to the broad maximum at T ~ 1 there is only one (or two) low-T maximum. 
Unfortunately, the finite-size effects are rather large for 0 < Ai < 1 and the temperature 
To is To ~ 0.01 and Tq ~ 0.1 for Ai = 0.5 and Ai = 0.05. This complicates considerably 
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the determination of the behavior in the thermodynamic limit. All the same, we believe 
that multi-peak structure for C'(T) or many stair-steps one for s(T) does not appear for 
0 < Ai < 1 and the model with Ai 1 is a generic one in the range 0 < Ai < 1. 

Thus, we suggest that the value Ai = 1 divides the regions of the parameter Ai with 
the qualitatively different behavior of the low-temperature thermodynamics. The possible 
argument in favor of this assumption is the fact that this model describes the transition line 
between the ground state phases of different types for Ai > 1 and Ai < 1 (see Figj2]). 


V. SUMMARY 

We have studied the ground state and the low-temperature thermodynamics of the delta- 
chain with anisotropic F and AF interactions. At dehnite relations between values of these 
interactions the lowest branch of the one-magnon states is dispersionless which means the 
existence of exact localized eigenfunctions. If the energy of the localized states is zero 
(lowest energy) then they are the ground states. In this case the model depends on a single 
parameter which can be chosen as the anisotropy of basal-apical spin interaction Ai. When 
Ai = cxD the model reduces to the Ising model with equal F and AF interactions while for 
Ai = 0 it is the isotropic ferromagnet on the basal chain and independent spins on the 
apical sites. Remarkably, the ground state degeneracy is the same for all 0 < Ai < oo 
(excluding some special values of Ai). The degeneracy is macroscopic and leads to hnite 
residual entropy at T = 0. In the limiting cases Ai = 0 and Ai = oo the ground state 
degeneracy is even larger and for hnite Ai it is partially lifted. 

For Ai 3> 1 the excitation spectrum has multi-scale structure and is rank-ordered in 
powers of small parameter Aj"^. The number of sections of the spectrum is equal to the 
number of triangles of the delta-chain and the energy of the levels in the m-th section 
E ~ /\- 2 m +2 (jYi = \ _ N/2). The origin of such exponentially low energy levels is the fact 

that the m-magnon bound complex in this system has the energy Em ~ A^ Each m-th 

section of the spectrum is responsible for the appearance of m-th peak in the dependence 
C{T). Thus, the number of the peaks in the specihc heat grows with the length of the chain. 
Similarly, such multi-scale structure of the spectrum determines the characteristic features 
of the dependence of s{T) and x(T). Numerical calculations by the ED of hnite delta-chains 
for smaller values of Ai (up to the isotropic point Ai = 1) show that the behavior of the 


24 


thermodynamic quantities is qualitatively similar to that for Ai S> 1. 

For Ai 1 the spectrum has two-scale structure in contrast with the case Ai S> 1. It 
leads to only one low-temperature maximum in C{T). According to the numerical calcula¬ 
tions such feature in the behavior of C{T) survives in the region Ai < 1. Thus, the isotropic 
point Ai = 1 separates the region of Ai on two parts with the qualitatively different behav¬ 
ior of low-temperature thermodynamics. We note that this conclusion is based on numerical 
calculations and needs more rigorous analysis. 

The main points of our study can be extended to delta-chains with arbitrary spin values 
because the condition on exchange interactions for existence of localized magnons does not 
depend on the value of spin s. 
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